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ABSTRACT 

We analyze subhalos in the Via Lactea II (VL2) cosmological simulation to look for corre- 
lations among their infall times and z = dynamical properties. We find that the present- 
day orbital energy is tightly correlated with the time at which subhalos last crossed into the 
virial radius. This energy-infall correlation provides a means to infer infall times for Milky 
Way satellite galaxies. Assuming that the Milky Way's assembly can be modeled by VL2, we 
show that the infall times of some satellites are well constrained given only their Galactocen- 
tric positions and line-of-sight velocities. The constraints sharpen for satellites with proper 
motion measurements. We find that Carina, Ursa Minor, and Sculptor were all accreted early, 
more than 8 Gyr ago. Five other dwarfs, including Sextans and Segue 1, are also probable 
early accreters, though with larger uncertainties. On the other extreme, Leo T is just falling 
into the Milky Way for the first time while Leo I fell in ~ 2 Gyr ago and is now climbing out 
of the Milky Way's potential after its first perigalacticon. The energies of several other dwarfs, 
including Fornax and Hercules, point to intermediate infall times, 2 — 8 Gyr ago. We compare 
our infall time estimates to published star formation histories and find hints of a dichotomy 
between ultrafaint and classical dwarfs. The classical dwarfs appear to have quenched star 
formation after infall but the ultrafaint dwarfs tend to be quenched long before infall, at least 
for the cases in which our uncertainties allow us to discern differences. Our analysis suggests 
that the Large Magellanic Cloud crossed inside the Milky Way virial radius recently, within 
the last ~ 4 billion years. 
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1 INTRODUCTION 

The Milky Way is a unique laboratory for understanding the lives 
of dwarf galaxies (L < 1O 8 M0). Dwarf spheroidal galaxies, in 
particular, stand out among galaxies because of their high dark 
matter content, lack of gas, and lack of recent star formation. Like 
larger galaxies dDressleJl98 0; Butche r"& OernleJl984l : lGoto et alj 
2003), dwarf galaxies appear to have a "morphology-density" rela- 
tion, with dwarf spheroidal galaxies preferentially crowding around 
norma l galaxies (or with in groups) instead of the field jMated 
Il998al : IWeisz et all 1201 lh . All of the galaxies within the Milky 
Way's dark-matter halo except the two Magellanic Clouds are 
dwarf spheroidals. 

In galaxy groups and clusters, the morphology-density rela- 
tion is associated with a color-density relation or a star-form ation- 
density relation jKauffmann et alj|2004 iBlanton et alj|2005t) . The 
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origin of these relations is thought to result from the quench- 
ing of star forma tion once galaxies become satellites in larger 
systems (see, e.g. iBerrier et alj|2009l) . Once inside the virial ra- 
dius of a larger host, star formation in the satellites may be 
quench ed either because they stop accre t ing fresh gas ("strangu- 
lation " ; ILarson. Tinslev & Caldwell 1980blBekki. Couch & Shioval 
l2002h or b ecause their cool g as is stripped away ("ram-pressure 
stripping"; iGunn & Gottll 19721) due to interactions with the host's 
gas halo. High-speed encounters with other satellite galaxies or 
the host itself may similarly affect mo rphologies and star forma- 
tion ("harassment"; iMoore et al. 1996). These processes are also 
relevant for dwarf galaxies around the Milky Way. The Milky 
Way is likely surrounded by a hot gas halo of its own, which 
can ai d in quenching star formation once galaxies fall within its 
reach jMaller & Bullockll2004l: iFang et alJl2006l:[Peek et alj|2007l : 
iKaufmann et al.ll2008l : iGrcevich & Putmanl 120091. l2010h . In addi- 
tion, dw arfs can experience dynamical transformatio n due to tidal 
stirring fcokas et al.l l20ld ; iKazantzidis et alj 1201 if) , though the 
smallest gal axies may very well be born puffy w ell before enter- 
ing the halo jKaufmann. Wheeler & Bullock 20071). 
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One way in which dwarf satellites are different than larger 
satellite galaxies in groups and clusters is that their low mass 
makes them inherently fragile, thus more susceptible to quench- 
ing processes that would not otherwise affect ~ L* galax- 
ies. Specifically, it is possible that some of the known dwarf 
spheroidals were quenched prior to infall into the Milky Way. 
During and after reionization, photons pour into the intergalac- 
tic medium, heating and pressurizing the gas so much that it is 
unable to collapse onto dark-matter halos with circular velocity 
smaller t han Vmax ~ 30 km s~ (the exact value of which is 
dispu t ed: iThoul & Weinberdll996l:lBullock, Kravtsov & Weinberg! 
200d : IBenson et al.l |2002|; IBenson & Madaul |2003l : iDiikstra et alj 
20041 ; lOkamoto. Gap & Theunsl l2008). Reionization photons may 
also ph otoevaporate gas alread y present in halos before reion- 
ization ( Barkan a & Loebl 1 1999h . Dwarf-galaxy dark-matter halos 
have small escape velocities; therefore, stellar winds or super- 
novae may permanent ly blow gas out of these small galaxies 
dGovernato et al.l 2010). The newly discover ed population of ul- 
trafaint dwarf spheroidals in the Mi l ky W ay Jwillman et al] |2005 ; 
iBelokurov etafl 120071 : iKirbv et al.l 120081) have overwhelmingly 
old stellar populations and are often speculated to be "fossils of 
reioni zation" — galaxies that only form stars prior to reionization 
(e.g.. iMartin. de Jong & Rb3l2008l ; |Madau 2009; iBovill & Ricottil 
l2010al lbl). However, it remains unclear whether the ultrafaints are 
old because they are true fossils of reionization or simply because 
they fell into the Milky Way at an early epoch. 

One way we can hope to discriminate quenching scenarios is 
to determine when each galaxy became a satellite and then com- 
pare this to an inferred star formation history. The Milky Way is a 
unique laboratory for answering these questions not only because it 
is currently the only place where we can find ultrafaint galaxies but 
because of the availability of exquisite photometric and kinematic 
data for virtually all of the satellites. We have three-dimensional 
configuration-space positions and line-of-sight velocities for ev- 
ery dwarf satellite (the dwarfs considered here are listed in Table 
1). Moreover, for a subset of the classical dwarf galaxies, we also 
have proper motions. Thus, we have four- or six-dimensional phase 
space positions for the Milky Way dwarf satellite galaxies. These 
kinematic data allow us to estimate the infall times of the dwarfs. 

Previous work has focused on constraining the dwarf infall 
times by evolving satellite orbits back in time based on those ob- 
served phase-space coordinates today or by tracing sp ecific satel- 
lite orbits forward in time in lar g e N-body simulations dBesla et alj 
20071: iLux, Read & LakdbOld: lAngus, Diaferio & Kroupj I20T1I ; 
Bovlan-Kolchin. Besla & Hernquistl20Hl) . The problem with these 
approaches is that they are sensitive to Poisson noise — specific 
things like the choice of the triaxiality of the Milky Way and its 
evolution through time, satellite interactions in the simulated Milky 
Way, treatment of dynamical friction and tidal stripping of the satel- 
lites, all cause large uncertainties for the infall times of individual 
orbits. 

Instead, we adopt a simpler, more statistical approach to deter- 
mining the infall times of the Milky Way dwarf satellites (including 
the dwarf irregular Magellanic Clouds). In particular, we focus on 
using a simulation of a Milky- Way-type halo to determine an infall- 
time probability distribution function (PDF) for each dwarf based 
on simulated subhalos with similar present-day phase-space coor- 
dinates. In Sec. [2] we describe the properties of the simulation that 
are relevant to this work. We use the simulation to show, in Sec. [3] 
that there is a strong correlation between the infall time of subhalos 
that might host dwarf galaxies to their binding energy today to the 
host dark-matter halo. We call this correlation the "energy-infall" 



relation. Since the simulated halo is similar to the halo that hosts 
the Milky Way, we make the ansatz that the energy-infall relation 
of Sec. [3] can be applied to Milky Way dwarf galaxies. In Sec. [4] 
we create PDFs of the dwarf infall times based on the subhalos that 
have galactocentric positions, line-of-sight velocities, and proper 
motions (if measured) within the measurement error bars of ob- 
served dwarf galaxies. We can constrain the infall times using the 
energy-infall relation because the kinematic data yield estimates 
of the binding energy (or upper limits thereof if no proper motion 
measurements exist). In Sec. [5] we show how the infall-time PDFs 
correspond to existing determinations of the star-formation histo- 
ries of the Milky Way's dwarfs. We discuss what the implications 
are for quenching mechanisms for dwarf spheroidals and trends be- 
tween galaxy properties and environment. In addition, we highlight 
further applications of the energy-infall relations in the study of 
galaxy evolution. We summarize our conclusions in Sec. [6] 



2 METHODS 

We use the Via Lactea II (VL 2) simulation for our analysis of 
satellite orbits and infall times (|Diemand. Kuhlen & Madad l 20071 : 
iDiemand et alj2008l : lKuhlenll2010l) . The VL2 simulation is a high- 
resolution lambda cold dark matter (ACDM) cosmological simula- 
tion that focuses on a dark matter halo of approximately the same 
size as the one that hosts the Milky Way, with a maximum circular 
velocity V max = 201 km/s at z = 0. 

The cosmology assumed in this simulation is taken from the 
flat- universe six-param eter analysis of the WMAP three-year data 
set dSpergel et al]|2007h : fi m = 0.238, fl A = 0.762, h = 0.73, 
n s = 0.951, as = 0.74. The resolution of VL2 is high enough, 
with particle masses of 4.1 x lO 3 Af0 each, to resolve thousands 
of subhalos bound to the main host. We have selected a sam- 
ple of ~ 2000 bound subhalos with maximum circular velocities 
Vmax > 5 km/s at z = and studied their dynamics at redshift 
z = 0. Hal os in VL2 are found through the 6DFOF halo finder de- 
scribed in lDiemand. Kuhlen & Madaul d2006h ■ The Vmax threshold 
guarantees that we do not miss any subhalos that might plausibly 
be hosting dwarf galaxies in the Milky Way today. 

Throughout this work, we define the main-halo mass A'hoo {z) 
and virial radius R2oo(z) in terms of an overdensity A = 
200/5 m (z), where p m {z) is the homogeneous matter density in 
the Universe, and M200 = 47r Ai?2oo /3- With this definition, the 
z = virial mass of the VL2 host is M 20 o = 1.9 x 10 12 M S today, 
with a virial radius of just over 400 kpc. We define the center of the 
halo to be the center of mass of bound particles within the virial 
radius. 

The infall time of a subhalo into the main host depends some- 
what on definition. We identify t in f a u as the look-back time since 
the subhalo last crossed inward through the virial radius of the main 
halo R,2qq{z). We explored an alternative definition where ti n f a ii 
is identified with first crossing through the virial radius but found 
that the resulting kinematic correlations were more scattered (but 
still present), so we adopted the last-crossing definition. 

We define the binding energy of a subhalo to the host halo by 

£ = -<t>(r)-±V\ (1) 
where the gravitational potential 

is defined such that 4>(Ro) = for a reference radial position Rq. 
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Figure 1. Binding energy vs. infall time for the selected sample of VL2 
subhalos at z = 0. Colors indicate galactocentric distance. Notice how the 
least bound subhalos are the ones accreted most recently and also the only 
ones with large galactocentric distances. 



The enclosed mass at galactocentric distance r is M(< r) and the 
subhalo velocity with respect to the halo center is V. For this work, 
we define Ro = 1 Mpc in physical units so that the gravitational 
potential has a fixed zero point across cosmic time. By contrast, if 
we had chosen Ro = R200, the energy of a particle would change 
throughout time even if its orbit were fixed and the density profile 
were constant in time since R200 grows with time. 



3 THE ENERGY-INFALL CORRELATION 

We checked for correlations among subhalo infall times and many 
different subhalo orbital properties (orbit circularity, angular mo- 
mentum, binding energy, radial velocity, current position, etc.). 
Many of these properties showed no strong correlation with infall 
time, but there is a clear correlation with binding energy. A key re- 
sult of this paper is that the current binding energy of a subhalo £ 
is a simple, clean 2 = indicator of the subhalo's infall time. 

The energy-infall correlation is demonstrated in Fig. [T] in 
which we also color-code the subhalos by their current radial posi- 
tion with respect to the halo center. We find that subhalos accreted 
early are more tightly bound to the halo with relatively little scatter, 
and that late-infall subhalos have lower binding energy (albeit with 
more scatter). Another way of looking at this correlation is that we 
find that subhalos that are currently deep within the potential well 
of the halo preferentially have higher binding energy and earlier 
infall times (greater tinfall) than subhalos that orbit farther out. 
While this is generally expected from hierarchical structure forma- 
tion theory, the quantification of this relation we find in this work 
proves highly useful, and we return to this point in later sections. 

We can trace this correlation directly to the energy of the sub- 
halos at their infall epochs. In Fig. [2] we show the ratio of the bind- 
ing energy of subhalos today to the binding energy at ti n f a ii, as a 
function of ti n f a ii- Recall that £ is larger for more tightly bound 
subhalos, so that the £i n faii/£today < 1 for halos that become 
more bound after infall and £i n faii/£today > I if the subhalo 
becomes less tightly bound towards z = 0. This figure demon- 



strates that while there is some scatter about £i n faii/£toda y = I, 
the average value of this ratio as a function of infall time is nearly 
f and does not change appreciably with infall time. Thus, the sub- 
halo on average conserve their energies at infall, with Emfall ~ 
—4>(R'2oo aU ) ~ 5(^200 , ! ) 2 . Because the binding energy at infall 
is linked to the virial properties of the host halo at that time, the 
correlation of binding energy with infall time arises from the mass- 
assembly history of the halo. 

Why are subhalo energies, on average, conserved throughout 
cosmic time? For the moment, we only consider changes of en- 
ergy due to interactions between subhalos and the host. Three-body 
interactions including another subhalo may also serve to increase 
subhalo energies, or reduce binding ener gies ( Sal es et a n 12007I : 
iD'Onghia & Lakel2008l ; lLudlow et ail2009n . but we consider only 
those processes that cause subhalos to become more bound. The en- 
ergy may change if dynamical friction is significant or if the grav- 
itational potential evolves. Let us consider the case of dynamical 
friction. Dynamical friction is only likely to change the energies of 
the most massive subhalos, since the merger timescale 



oc {M 2 oa/M sat ) 1/A t dyn , 



(3) 



where tdyn ls the typical dynamical time in the host 
alaxy (mass M200), and the s atellite mass is M sa t 
Bovlan-Kolchin. Ma & Ouataertl f2008). The merger timescale is 
linked to the timescale over which the energy changes, since it 
takes a time T merge for the satellite to go from an initial binding 
energy £ in f a ii to £ to da y = —0(0). For a Milky Way-mass halo, 
the dynamical time is roughly 1 Gyr, and so only the subhalos 
that have M sa t/M2oo > 0.1 will have merged with the host 
halo in a Hubble time. Thus for most z — subhalos, dynamical 
friction will have only a modest effect on the binding energy, in 
the direction of making them increasingly more bound. Those 
subhalos for which dynamical friction is important in changing 
the energy significantly are also those that are most likely to have 
already merged with the host or have been tidally shredded, and 
are thus not part of the surviving z = subhalo population. 

Let us consider the second case of changes to the halo poten- 
tial. Dark-matter halos typically have a "fast" growth stage, within 
which the matter within the scale radius of the density profile to- 
day is rapidly acquired, and a "slow" growth stage, after which the 
halo g rows constantly a nd without major mergers dWechsler et alj 
2002; Zhao et al, 2009). These two regimes are artifacts of the 
shape of the linear density perturbation, with the fast growth regime 
linked to regions in which the density perturbation 8 ~ const, 
and the slow growth regime to regions in which the 5 is a sharply 
falling function of distance from the center of the perturbation 
dDalal. Lithwick & Kuhlerfeoid) . For most of the VL2 halo's his- 
tory, it is in the slow growth period, meaning that the density pro- 
file near and outside the virial radius is a sharply falling function 
of distance. If we approximate the local density profile as a power 
law, then p oc r~ a , the enclosed mass M(< r) oc r 3 ~ a and the 
gravitational potential 



R " GM(< r) 



Rl~ a ~r 2 - a . 



(4) 

(5) 
(6) 



Thus if the density profile falls off as a > 2, which is certainly the 
case in the slow-growth phase, then 4>{r) oc r 2 ~ a . Since halos form 
from the inside out, the enclosed mass at a fixed radius within the 
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Figure 2. Ratio of the binding energy at infall to the binding energy today 
as a function of subhalo infall times. On average the subhalos in the sample 
are slightly more bound today than when they fell into the VL2 halo, but 
clearly the binding energy of subhalos today tracks the binding energy at 
infall regardless of the infall time. 



virial radius does not change much with time teook et alJ|20lTh . 
As such, we do not expect the potential to change much with time 
at a fixed radius. 

In summary, we find a strong correlation between the energy 
of subhalos at z = in the VL2 simulation and their infall times. 
This is because the subhalos surviving to z = on average con- 
serve their binding energies since their infall epoch, and all halos 
falling in at the same time must have roughly similar energies, of 
order the energy of a circular orbit at the virial radius. 



4 INFALL TIMES OF MILKY WAY SATELLITES 

Under the assumption that Milky Way satellite galaxies are hosted 
by dark matter subhalos similar to those predicted by the VL2 simu- 
lation, and that the Milky Way itself and its mass-assembly history 
is something like the VL2 halo, we can estimate the infall times 
of Milky Way satellites based on their energies, which rely funda- 
mentally on kinematic measurements. The goal of this section is to 
show how the measured kinematics of the Milky Way dwarfs define 
a probability distribution function (PDF) for the infall time of each 
dwarf. Ideally, we would use positions, line-of-sight velocities, and 
proper motions for each dwarf so that we could reconstruct each 
dwarf's orbital energy assuming that the gravitational potential of 
the Milky Way matches that of VL2. We would not expect a delta- 
function like peak in the PDF for t in f a u even in this idealized case 
due to the scatter in the energy-infall relation, but we would expect 
a well-defined peak in the PDF (see Fig. []}. From Fig. Q] we see 
that the position information should reduce the width of the PDF, 
especially for recently accreted subhalos. 

Unfortunately, for most Milky Way dwarf galaxies, especially 
the subset of ultrafaint dwarfs, we only have position measurements 
and line-of-sight velocities. Thus, we begin by exploring how well 
one may constrain the infall time with only these four-dimensional 
(position vector and line-of-sight velocity) data. 



4.1 Constraints from radial velocities and distances alone 

In Fig. [3] we present our subhalo population in the space of Galac- 
tocentric distance r and radial velocity V r . The subhalo points are 
color coded by infall time, as indicated by the legend on the right. 
The outer envelope in r — V r space is dominated by recently ac- 
creted subhalos; as one moves to smaller \ V r \ and especially as one 
moves to smaller r, the subhalos are accreted further back in time. 
The red points (ti n f a ii ~ 1 Gyr) with negative velocities and large 
radii correspond to systems that are just falling in for the first time 
and the yellow points (t in f a u ~ 3 Gyr) are systems that are just 
coming back out after their first pericenter passage. The purple and 
blue points correspond to early accretions (tinfall ~ 8 — 10 Gyr) 
and naturally cluster at the small radii and fairly low speeds indica- 
tive of higher binding energy. 

The right-hand panel of Fig. [3] is identical to the left except 
that we have overlayed the Galactocentric distances and line-of- 
sight velocities of Milky Way dwarfs. The photometric and kine- 
matic data on the Milky Way dwarf galaxies we consider are pre- 
sented in Table 1. We identify the measured line-of-sight veloc- 
ity for the dwarfs with the radial velocity since the Sun is deep 
in the potential well of the Milky Way and the dwarfs are much 
farther out. We can estimate the infall times of dwarfs by compar- 
ing their positions on this diagram to the infall times of the VL2 
subhalos that occupy similar positions. For example, we can dis- 
cern that the Leo T and Leo I dwarf spheroidals were accreted 
recently, with t in f a u < 1 Gyr and Unfaii ~ 2 Gyr respectively. 
We note that an infall ti me of ~ 2 Gyr for Leo I is consi stent 
with the model reported bv lMateo. Olszewski & W alker ( 20081) that 
gives the last pericenter crossing for Leo I as ~ 1 Gyr. Conversely, 
Segue 1 appears to be an old companion of the Milky Way with 
tinfaii ~ 8 — 10 Gyr. The infall times for other MW satellites are 
not as easily discerned by eye in this diagram, but many still pro- 
vide useful constraints on the infall time even in the absence of 3d 
velocity data, as we will now discuss. 

The solid black histograms in Fig. [4] show infall-time PDFs 
for most of the known MW satellites. We constructed these based 
on each galaxy's r and V r values compared to those of the VL2 
subhalos. The top three rows show the classical dwarf galaxies and 
the bottom four show the more newly discovered ultrafaint dwarf 
population. The red histograms include proper motion information 
and will be discussed in Sec. 14.21 

The solid black PDFs for each dwarf d are constructed by in- 
cluding the infall times for all VL2 subhalos with radii r that obey 
\r — r d \ < c ay and radial velocities V r that obey \V r — V*\ < 
c av r ■ Here oy and eryy are the observational errors on radial ve- 
locity and galactocentric distance for the appropriate galaxy d. The 
variable c > 1 is a proportionality constant set by the requirement 
that there are at least 20 VL2 subhaloes in the subsample of each 
galaxy. The positions, radial velocities, and uncertainties for the 
Milky Way satellites are summarized in Table |5.l31 

Based on their radial positions and velocities alone several of 
the dwarfs have reasonably well constrained infall times. As previ- 
ously discussed, Leo I and Leo T have been quite recently accreted. 
Draco, Sextans, and Bootes likely fell into the Milky Way well be- 
fore z = 1 (8 Gyr), although the precise redshift beyond z = 1 is 
not clear. The infall times for other dwarfs are ambiguous. In par- 
ticular, Canis Venatici II and Willman I could have been accreted 
anytime in the past 10 Gyr. 
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Figure 3. Infall times in color as a function of the radial velocity and the galactocetric distance of the VL2 sample of subhalos. The left panel shows only 
the VL2 subhalos, the right panel shows a sample of Milky Way dwarf galaxies on top of the VL2 subhalos. It is evident that with just the radial velocity and 
galactocentric distance information coupled with the VL2 predictions we can set constraints on the infall times of Milky Way dwarfs from their position on 
this space. It is specially clear that LeoT and Leol seem to have been accreted to the Milky Way halo recently as opposed to Segue 1, which appears to be 
accreted early. 



4.2 Including proper motions 

Six MW dwarfs have published proper motion estimates from 
the Hubble Space Telescope (HST) with reasonably small errors: 
Ursa Minor, Carina, Sculptor, Fornax, the Small Magellanic Cloud 
(SMC), and the Large Magellanic Cloud (LMC; see Table 1, 
which provides associated references in each case). Leo II also 
has a published pro per motion, but its relative error is fairly large 
(Lepin e~et alj|20lih . For the galaxies with small errors we should 
be able to to place better constraints on their accretion times as we 
may estimate directly the binding energy and not just place an up- 
per limit on it. 

To demonstrate the improvement to the infall time constraints 
coming from proper motions we examine a subsample of VL2 sub- 
halos with similar radial velocities and galctocentric distances to 
each of the MW dwarfs under study. Again, the subhalo sample for 
each Milky Way satellite is selected such that there are at least 20 
subhalos in each sample according to the position and radial veloc- 
ity constraints discussed in the previous section. Figure|5]shows the 
tangential velocities and infall times of those subhalos associated 
with each of those satellite galaxies for which proper motions are 
measured with high precision. The grey band shows the la mea- 
surement of the proper motion for each dwarf. It is clear from Fig. 
[5]that the extra information added from proper motions can sharpen 
the infall-time PDFs. 

The red histograms in the upper panels of Fig. [4] demonstrate 
that the extra information does indeed sharpen the infall-time con- 
straint relative to the constraint without proper motions. For Ca- 
rina, Ursa Minor, and Sculptor, their relatively small tangential 
velocities disfavor the tail of recent infall-times that would be al- 
lowed by their radial velocities alone (compare red to black solid 



histograms), pushing their inferred infall times towards earlier ac- 
cretion (tmfaii > 8 Gyr). For the SMC, on the other hand, its 
relatively high proper motion favors later infall solutions, giving 
tinfaii < 4 — 9 Gyr. Interesting, these late infall solutions for 
the SMC were otherwise disfavored based on its radial velocity 
alone. A similar story follows for the LMC, which must be ac- 
creted recently according to our comparison with VL2 subhalos 
(tinfaii i$ 4 Gyr). The measured tangential speed for Fornax is 
low enough to disfavor the possibility of very recent accretion, giv- 
ing ti„faii ~ 5 - 8 Gyr. 

Finally, the proper motions reported for Leo II provide little 
additional constraint on its infall time, though reassuringly Leo II 
proper motions are consistent with those expected for subhalos of 
the appropriate distance and radial velocity. It should be empha- 
sized that the same agreement applies to the SMC and the LMC. 
We do find VL2 subhalos in our sample that can be matched with 
their speeds and positions. 

Table 1 provides a summary of our results on infall times in 
the right-most column. 

4.3 Future observations 

Since a number of the dwarfs without proper-motion measurements 
at the present, especially the ultrafaint population, have broad 
infall-time PDFs, we would like to know how well one would need 
to measure the proper motions to get unambiguous infall times. 
In Fig. [6] we show scatter plots for the tangential velocity as a 
function of infall time for the VL2 subhalo samples associated 
with Ursa Major I and Willman I. The nature of t he latter object 
is und er debate; it may be a disrupting star cluster Iwillman et alj 
However, the distribution of tangential velocities as a func- 
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Figure 4. Infall time histograms for each dwarf constructed by matching kinematic properties of subhalos in VL2 to measured properties in the dwarfs (see 
text for details). The solid black lines include only radial velocity and position information for the dwarfs. The red histograms add the additional constraint 
associated with proper motions when they are available. When proper motions errors are relatively small (top two rows) this additional information tightens 
the inferred range of infall times significantly. 



tion of infall time for associated subhalos is illustrative. For both of 
these objects, if the tangential velocity is small (< 200 km s _1 ), a 
proper-motion measurement with an associated uncertainty in the 
tangential velocity of ~ 50 km s" 1 (typical of the uncertainty in 
the tangential velocity of the classical dwarfs) would be sufficient 
to bracket the infall time to within 2 Gyr. The scatter in tinfaii for 
fixed tangential velocity is higher for large tangential velocities, but 
even an uncertainty of 50 km s _1 would be sufficient to tell if these 
objects fell in early (t in f a u > 8 Gyr) or late. 



5 DISCUSSION 



Here, we discuss the energy-infall relation and infall-time PDFs 
in several contexts. First, we compare the infall-time PDFs for in- 
dividual dwarf galaxies with their star-formation histories. Second, 
we speculate as to what kinds of theoretical and observational work 
would be required to produce more robust infall-time PDFs. Third, 
we speculate on some possible applications of the energy-infall re- 
lation to tidal streams. 
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Figure 5. Tangential velocity as a function of infall time for subsamples of subhalos with similar radial velocities and galactocentric distances to those of the 
given dwarf galaxies. The subsample selection criterion is the same as in Fig. [4] The la uncertainties in the proper motions are represented by the shaded 
regions. The addition of proper motion constraints provides a better estimate of the infall time than radial velocity alone. 



5.1 Dwarf galaxy infall times and their star-formation 
histories 

One of the main reasons for determining the infall times of the 
Milky Way dwarf galaxies is to use that time with respect to 
the star-formation history in the dwarfs as a diagnostic for star- 
formation quenching mechanisms in these small galaxies. Al- 
though the infall-time PDFs are based on a single simulation of 
a Milky-Way-like halo, we show that the infall-time PDFs have the 
power to show a number of interesting trends regarding quenching 
in small galaxies. For this discussion, it will be helpful to the reader 
to have Fig.|4]in view or to refer to the last column in Table 1 , where 
we have compiled our infall time estimates. 



5.1.1 Classical dwarfs 

Among the classical dwarf spheroidal galaxies, we find that Ca- 
rina, Draco, Ursa Minor (UMi), and Sculptor likely fell into 
the Milky Way early (ti n f a ii > 8 Gyr). Of those four, UMi 
and Sculptor have old stellar populations, with no evidence of 
stars younger than t ~ 10 Gyr in UMi and t ~ 7 Gyr in 



Sculptor (Hurlev-Keller, Mateo & Grebel 1999; Monkiewicz et al. 


1999]: Mighell & Burkej 1999: Wyse et alj 


19991: ICarrera et all 


2003; Dotohirj|2002l: iTolstov et al.N2004ICc 


>hen & Huana l20ld: 
jminantly old stellar 


de Boer et al. 2011). Draco likelv has a predi 



most a uthors suggest most stars are > 10 Gyr old jGrillmair et al 



19981 : lAparicio. Carrera & Martinez-Delgadol l200ll ; lOrban et al 



2008). The Carina dwarf is somewhat different. It likely formed 



most of its stars in bursts at 3, 7 and 1 1 Gyr lookback time, with 
~ 50% of the stars 
1996 



7 Gyr ago ( Smecker-Hane et al. 
Hurley-Keller. Mateo & Nemec 1998); iMonelli etal 



1994 



2003 



Rizzi et alj |2003l) — after the time we estimate that it fell into the 
Milky Way. 

We find that three classical dwarf spheroids should have fallen 
in significantly later. Leo I is the most recent at ~ 2 Gyr ago, 
followed by Leo II at ~ 2 — 6 Gyr and Fornax at ~ 5 — 8 
Gyr look back time. Interestingly, all three of our late infall can- 
didates demonstrate active star formation at later times than do 
our early-infall candidates. Leo I & II both have long epochs 
of star formation stretching from > 1 Gyr ago to ~ 2 Gyr 



ng 

ago dMighell~&R ich 1996; Caputo et al. 1999; Gallart et alj[l999l: 
| Dolphiiil2002l;|Bosler. Smecker-Hane & Stetsoj2007l;|Orban et alj 
l2008l : lGullieuszik et alj|2009l : lLanfranchi & Matteuccill2010h . The 
end of star formation correlates well with the predicted infall time 
for Leo I and is consistent with the fairly broad infall range for 
Leo II. Fornax had a major burst in star formation 3 — 4 Gyr ago 
( Stetson, Hesser & Smecker-Hane 1998; Saviane, Held & Bertelli 



2000; Coleman & de Jong 2008; 



Held et alj I2TJ10I ; Letarte et al. 



2010), which corresponds to a time just after our estimated infall. 



population as well, and although the exact age of this is debated, 



Both the LMC and SMC are still active ly forming stars 
dPe Marchi et alj l201ll ; lHarris & Zaritskvl I2009T) and are unique 
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Figure 6. Tangential velocity as a function of ti n f a u for VL2 subhalos that have similar positions and radial velocities of Ursa Major I (left) and Willman I 
(right). Measured proper motions with an uncertainty of order ~ 50 km s — 1 can bracket infall times to within ~ 2 Gyr in most cases. 



among the classical Milky Way satellites in containing neutral gas 
today dGrcevich & Putmanl [20091 l2010h . According to our VL2 
subhalo analysis, the high 3D speed and position of the LMC de- 
mand that its accretion was more recent than ~ 4 Gyr. Similarly, 
while the radial velocity and position of the SMC alone would seem 
to favor an early infall (see Fig.[4jl its relatively high proper motion 
pushes us towards recent accretion < 4 — 9 Gyr albeit with large 
uncertainties. Though we have not demanded that subhalos be in- 
teracting or massive in order to assocate them with the Clouds, it is 
at least encouraging that our analysis is consistent with the idea that 
they were accreted fairly recently and a t approximately the same 
time jBesla et alJl2007l ; ISales et alj201 ll) . 

Fig.|7]provides an illustration of the diversity of star-formation 
histories relative to the infall times in the classical dwarfs exempli- 
fied by a few specific cases (with names indicated in each panel). 
UMi, Sculptor, and Leo II are examples of dwarfs that appear 
to be quenched just as they fell into the Milky Way halo. While 
star formation in Sculptor (upper middle) goes out with a whis- 
per, Leo I (lower middle) may or may not have experienced a 
burst of star formation after it fell in. The remaining dwarfs, il- 
lustrated by Carina and Fornax in the right panels, seem to expe- 
rience bursts of star formation immediately after infall. In galaxy 
groups and clusters, many galaxies experience bursts of star forma- 
tion triggered by the ho t host halo compressing gas in the satellites 
dGavazzi& J affe 198 5j; iBothun & Dres sier 1986 ; iBekki & Couchl 
20031: iMahaian Haines & Ravchaudhurvl l2010j: iRose et alj|2olol ; 
Abramson et alj|201 It ISantiago-Figueroa et alj 1201 11) . Something 



similar may be occuring with the dwarf galaxies in the Milky Way 
halo. 

At this point, it is not clear exactly what terminates the star for- 
mation or what the time scale for termination is because the infall- 
time PDFs are fairly broad, but it is useful to speculate in order to 
define a starting point. We suggest that the systems of the first cat- 
egory - the ones that that seem to have star formation quenched 



just upon accretion - are dominantly quen ched by quick-acting 
ram pressure stripping dGunn & Gotdll972h . The second class of 
objects, those with more gradual qu enching after accretion, are 
more likely affected by strangulation dLarson. Tinslev & Caldwelll 
1980). Finally, the bursting population may be so gas rich upon in- 
fall that they carry their material inward towards the Galactic center 
long enough to experience tidal shocking and associated bursting 
activity. These scenarios are clearly simplified, as all of these pro- 
cesses should operate to some extent and it is not clear why one 
effect should dominate over another from case to case. It may have 
to do with the details of the orbits or perhaps the arrangement of 
material within the infalling galaxy. Nevertheless our infall time 
estimates provide an initial point of comparison. 

5.1.2 Ultrafaint dwarfs 

Are the patterns any different for the ultrafaint dwarfs? One obvi- 
ous difference between the ultrafaint and classical dwarf popula- 



older, with the exception of Leo T dLee et al. 


120031: llrwin et alj 
Kuehn et alj2008l: 


2007l;lde Jong et al.l2008l|bllMartin et alj2008; 


Okamoto, Gao & Theuns 2008; Okamoto et al. 


20081: iNorriset alj 


2008; iGreco et alj|2008l; fevan- Weber et al.ll200 


8l; iLee et al.l 


2009|; 


Sandetal. 2009L l201fJ: ISimon et alj 20ld: Aden et all 


201 it 



Major II, Canis Venatici I and Leo IV which might have a small 
young (~ 2 Gyr) component dde Jong et alj l2008bl ; ISand et al] 
2010). Nevertheless, if the correlation between the end of star for- 
mation and the infall time characteristic of the classical dwarfs (and 
their quenching mechanisms) persisted for the ultrafaint dwarfs, we 
would expect the dwarfs to have disproportionally large Unfalh 

Indeed, we find that a higher percentage of the ultrafaint 
dwarfs we consider in this paper are consistent with an early in- 
fall (9/12; Sextans, Canes Venatici II Bootes I, Coma Berenices, 
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Figure 7. Qualitative depictions of star-formation histories (black solid lines) compared to our estimated infall times (cyan bands) for a sample of classical 
dwarf galaxies. The line heights in the Carina plot indicate the relative s t rength s of the star bursts. The like ly infall times are denoted by the shaded regions. 
The star-formation histories come from lHurlev-Keller. Mateo & Nemed jl998l) for Carina; iDolphii] (2002) for Ursa Minor, Sculptor, Leo I, and Leo II; and 
IColeman & de Jond f2003) for Fornax. 



Segue I, Ursa Major I, Ursa Major II, Willman I) compared to 
the classical dwarfs (5/9), where again we define early infall to be 
tinfaii > 8 Gyr. Yet the errors on the ultrafaint infall times are of- 
ten large enough that intermediate-age infall times are allowed for 
many of them. Only three ultrafaint dwarfs (Ursa Major II, Bootes 
I, and Segue 1) have kinematics and positions that actually demand 
early infall. 

Possibly more interesting are the ultrafaint dwarfs that were 
likely accreted somewhat later. Leo T is an outlier in the sense that 
it appears to be making its first passage into the Milky Way and it 
is the only Milky Way dwarf besides the Clouds to have a detection 
in neutral hydrogen dRvan- Weber et al.ll2008l ; iGrcevich & Putmanl 
120091. l20ld) . IGrcevich & Putmanl d2009l) find that Local Group 
dwarfs with galactocentric distances to either the Milky Way or 
M31 less than 270 kpc do not have detections in neutral hydro- 
gen while those at greater distances have a neutral gass mass 
> 10 5 Mq. Leo T, at a Galactocentric distance of ~ 410 kpc, is 
clearly a galaxy that has yet to be significantly influenced by the 
Milky Way. 

There are interesting differences among the other later-infall 
ultrafaint dwarfs. Canis Venatici I (ti n f a u ~ 2 — 6 Gyr) and Leo IV 
(~ 5 — 9 Gyr) both have predominantly old stellar po pulations but 
both show evidence fo r star formation ~ 2 Gyr ago jMartin et alj 
l2008l ; ISand et alj|20ld) . Unless the infall time lies at the wings of 



the infall-time PDFs for these objects, these galaxies were form- 
ing tiny numbers of stars (2% to 5% of the total) after they fell 
into the Milky Way. This is puzzling; a number of classical dwarfs 
continue to form significant numbers of stars up until ~ 1 Gyr af- 
ter infall (except for Carina and the Clouds, which keep bursting), 
and are continuing to form stars from early times until quenching 
inside the Milky Way halo. It is curious that Canis Venatici I and 
Leo IV both form a small number of stars after infall but that there 
was a long time gap between the formation of the ancient, dom- 
inant stellar population and the tiny young population. Infall oc- 
curred sometime in this gap. It means that although star formation 
was largely quenched before infall, these dwarfs held onto some 
gas that was prevented from forming molecular hydrogen until the 
dwarfs fell into the Milky Way halo. Perhaps additional pressure 
from the Milky Way hot halo was necessary to trigger the forma- 
tion of molecular hydrogen in those dwarfs. In any case, the star- 
formation histories of Canis Venatici I and Leo IV are very much 
unlike the star-formation histories of the more recently accreted 
classical dwarfs. 

Another interesting outlier is Hercules. T his dwarf galaxy 
shows no evidenc e for stars younger than 10 Gyr ( Sand et al. 2009; 
lAden et af]|201 lh . yet its infall time is predicted to be ti n f a ii ~ 
2 — 8 Gyr. This is the most significant offset we have found between 
probable infall time and the end of star formation other than the 
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Clouds. If our estimates of the infall time are accurate, it suggests 
that star-formation quenching must occur before the galaxy falls 
into the Milky Way, which is different than what we find for the 
classical dwarfs. This is different from all the other dwarf galaxies, 
for which we either cannot determine the order of star-formation 
quenching and infall or for which the end of star formation hap- 
pens after infall. 

This offset between the star-formation epoch and infall time 
of Hercules is even more interesting when one considers Leo IV. 
Leo IV and Hercules have quite similar Galactocentric distances 
(~ 140 — 150 kpc) and luminosities (~ lO 4 !/©), and both ap- 
pear to have fallen into the Milky Way since z = 1 (see Ta- 
ble 1). Neither has any detectable neutral gas. However, there is 
evidence for a small and yo ung stellar population i n Leo IV but 
Hercules is entirely ancient dSand et al"]|2009l l2010h . Why is this 
the case? One possibility is because the central dark matter den- 
sity of Hercules appears to be lower than al l the other Milky Way 
dwarfs JStrigari et al]|2008l : IXden et al.ll2009l) . which would place 
in among the smallest dark matter halos known to host a galaxy. 
The difference in the star-formation history of these two otherwise 
similar systems may hint that Hercules represents some sort of tran- 
sition in star formation or stochasticity in star formation in small 
objects. Better infall-time PDFs might shed light on the origin of 
the differences in the stellar populations given the similarity other- 
wise of these two dwarf galaxies. 

5. 1.3 Trends and Future Discovery Potential 

So far, we have found several fundamental differences in the 
quenching of classical and ultrafaint dwarfs and one possible simi- 
larity. In the cases in which we can distinguish the infall from star- 
formation epochs, the classical dwarfs appear to be quenched after 
infall but the ultrafaint dwarfs tend to be quenched for the most 
part long before infall (though uncertainties are large). This may 
suggest that the morphology-density relation or star-formation- 
density relation do not exist in the same way for dwarf galax- 
ies as they do for L, galaxies. While quenching in the classical 
dwarfs appears to occur within the Milky Way halo (and hence in- 
duces a star-formation-density relation), for the ultrafaint dwarfs 
for which the infall time can be cleanly separated from the star- 
formation epoch, the quenching appears to precede infall. This 
means that there is only a star-formation-density relation in that 
the reionization epoch is density-dependent, and various processes 
associate d with reionization can quench star formation in small 
galaxie s I Bulloc k, Kravtsov & We inberg 2000; B enson et alj2002j ; 
iBusha et alj20ld) . The one similarity is that it appears that a burst 
of star formation could be triggered in both classes of dwarfs after 
infall, but the degree to which that happens for both populations, 
and what that says about the specific quenching mechanisms after 
infall, are debatable. 

An interesting question to ask is what types of dwarf galax- 
ies we may expect to discover in next-generation deep wide-field 
surveys (e.g., LSSlQ), 

and what we may learn about star formation 
in these smallest of galaxies. Since the sample of bright (classical) 
Milky Way dwarf galaxies is likely complete, the question really 
revolves around the population of Milky Way ultrafaint galaxies, 
and how representative the current sample is of the population as a 
whole. One of the issues with the currently known sample of Milky 
Way dwarf galaxies is that we currently can only find ultrafaint 

1 http://www.lsst.org 



dwarf galaxies that are relatively nearby unless they are also rel- 
atively bright. From Fig. [3] we find that nearby galaxies dispro- 
portionately fell into the Milky Way early. Thus, we expect that 
next-generation surveys, which will be complete to the Milky Way 
virial radius for dwarfs with the same surface brightness as already- 
discovered dwarfs, could find a numbe r of faint dwarf galax ies that 
fell into the Milky Way more recently l lTollerud et alj |2008). 

There are already hints of interesting things in the current 
sample of dwarf galaxies, and it will be highly interesting to see 
if those hints become real trends when the catalog of Milky Way 
dwarfs is more complete. Since most of the ultrafaints already dis- 
covered are close, and tend to have early infall times, it is difficult 
to tell whether star-formation quenching preceded or followed in- 
fall. In other words, it is difficult to tell if the early-infall dwarfs are 
old because they fell into the Milky Way halo early and were then 
quenched, or if they were "born old" before they were accreted by 
the Milky Way. It may be easier to tell the sequence of events for 
more recently accreted dwarfs since the separation of time scales 
may be larger. If we find a number of more recent arrivals have 
only old stellar populations like Hercules, it will indicate either that 
feedback from those early epochs of star formation drove out all the 
gas in the galaxies or that the host halos were unable to accrete gas 
from the IGM. This clearly did not happen with Leo T, but this is a 
fairly bright galaxy compared to something like Segue 1 . 

If we discover new arrivals that are more similar to Canis Ve- 
natici I or Leo IV in that they have some young stars in addition 
to the predominantly old stellar population, we will learn that even 
the smallest galaxies may retain enough gas to late times to form 
some stars, even if formation of molecular hydrogen for star for- 
mation is suppressed. It would be surprising but highly interest- 
ing if we found ultrafaint dwarf galaxies that were not dominated 
by old stellar populations, since all of the currently known ultra- 
faint dwarfs have overwhelmingly old stellar populations. A num- 
ber of classical dwarfs are dominated by intermediate-age popula- 
tions (e.g., Fornax, Leo I, Leo II, Carina), and we should be able to 
leam at what mass-scale (either in stellar or dark-matter mass) these 
intermediate-age systems disappear if they do at all. One worry to 
keep in mind is that the lowest mass dark matter halos will tend to 
host dwarfs that a re the lowest surface brightness, and therefore the 
hardest to detect (Bul lock et alj|2010l) . Nevertheless, if we see any 
diversity in the star-formation histories for recently accreted ultra- 
faint dwarf galaxies, it will tell us something about the stochasticity 
of star formation in small galaxies and possibly about the physical 
mechanism(s) thereof. 

The questions we will likely be able to answer with a more 
complete catalog of Milky Way dwarf galaxies are the following: 
Are old galaxies old because they were quenched when they fell 
into the Milky Way or because they were unable to accrete or hang 
onto their gas prior to infall? How much is star formation quenched 
in small galaxies before they fall into the Milky Way? What prop- 
erty of the satellites is most correlated with quenching before in- 
fall, and what will that tell us about the physical mechanism(s) for 
quenching before the satellites fall into the Galaxy? How is star 
formation quenched once satellites are within the virial radius of 
the Milky Way? Do the satellites at larger radius have later epochs 
of star formation due to their infall time or because quenching is 
more effective at smaller galactocentric distances? How important 
is ram-pressure stripping for quenching star formation relative to 
other environmental factors? Can we learn something about the 
Milky Way's hot halo of gas today or its evolution in the past? How 
much does any of this depend on the mass of the dark-matter halos 
of the satellites? 
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Table 1: Observed and derived properties of the Milky Way dwarf satellite galaxies considered in this paper. 
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Note: Galaxies are grouped from top to bottom as pre-SDSS/classical MW dSphs followed by post-SDSS MW dSphs, with the Magellanic Clouds at last. 
Columns 5-7 show radial, tangential and spacial velocities in the galactic rest frame for those galaxies for which proper motions are known. For those 
galaxies with unknown proper motions column 5 shows the line of sight velocity in the Galactic Standard of Rest (GSR) frame yj os (gsr) . 
References: Except for Hercules and Leo T, values in column 3 (Luminosity) of the post-SDSS MW dSphs are from M artin, de Jong & Rix 12008). Values in 
columns 4 (Galactocentric Radi us) and 5 (Vj OB (gsr)) ar e derived and quoted respectively from the NASA/IPAC Extragalact ic Database. T he individual 
references are as follows: a) Pietrzvnski et al. (2009), b) D eriv ed from apparent mag nit udes liste d in Mateo 1 1998b), c) Bonanos et al. (2004), d) 
iBellazzini et al.|j200 4), e) Bell azzini. Gennari & Ferrarol l200$). fllColeman et al.U2007l). g)|Pietrzvhski et alU2008l). h)lLee et al.H2003IU)ICarrera et alj 
<2002h. illDall'Ora et alj J2006h. k ^ lMartin et all 120081). I l lGreco et all 12008I). m)|Belpkurov et al] )2007l) . n )ISand et all J2009h . o) lde Jong et all |20O83) ~d) 
lOkamoto. Gao & Theunsl 120081). q)IZucker et all <2006l). rtlWillman et alj J2005I). s)|Piatek et alj i2003l), t)|Piatek et all |2007l) , u) lLepine et aTffeOllI) , v) 
IPiatek et alj feOOfj) . wl lPiatek et alj <2005l) . x) lKallivavalil et alj j2006ir v) lKallivavalil. van der Marel & Alcockl 12006) . z) NASA/IPAC Extragalactic 

Database . 



5.2 Towards more robust and accurate infall-time PDFs 

Of course, the caveat to the discussion so far is that we have 
only examined the energy-infall relation and its consequences 
in one simulated dark-matter halo, and a halo simulated with 
(78 smaller by ~ 0.07 from the current preferred value. More- 
over, the actual mass of the Milky Way dark-matter halo is not 
known to better than a factor of two, so it is not clear if VL2 
is even mass-wise (let alone mass-assembly-wise) a good match 



ily - 

to the Milky Way ^Wilkinson & Evanj[T99a:lBattaglia et alj 



2006: 



Dehnen. McLaughlin & Sachania l l200^ T lXue et al 



2005 



2008 



Rcid et al. I l2009l; Watkins, Evans & An 2010). There is also some 
nontrivi al uncertainty in the mas s and distribution of stars in the 
Galaxy teinnev & Tremainell2008l) . The VL2 mass is on the high 
side, but certainly within range of what is expected for the Milky 
Way, given the existence of the Magellanic Clouds, which push 
expectations towards the ~ 2 x 10 12 Mq virial mass range 
( iBovlan-Kolchin. Besla & Hernquistll201ll) . 

There are several things that are likely to matter for getting 
the energy-infall relation right for the Milky Way specifically. The 
most important things may be both the dark and baryonic masses 
of the Milky Way, and the evolution of the baryonic component. 
Both the dark and baryonic masses affect the gravitational potential 
and hence the energy normalization in the energy-infall relation. In 
addition, the normalization of the subhalo mass function and dis- 



tribution in the halo depend on t he dark and baryonic masses of 
the Galaxy ( D'Onghia et al. 2010). The shape of the subhalo mass 
function depends on the baryons, too, since tidal shocking by the 
baryonic disk eliminates many of the lower-mass subhalos that oth- 
erwise would have bee n relatively un affected by dynamical fric- 
tion and tidal stripping to'Onghia etai]|2010|) . The differences in 
the radial distribution and mass function of satellites may shift the 
infall-time PDFs at a given position in phase space. 

Getting erg might matter as well. Structure forms earlier for 
higher erg, and subhalos are accreted earlier for the currently ac- 
cepted value of (78 rather than that used for VL2. This might shift 
the slope of the energy-infall relation. However, the gravitational 
potential of the host may also be different because of as, and hence 
the typical subhalo energy at infall may also be different. It will be 
important to see if the energy-infall relation depends on the under- 
lying cosmology. 

While it would be good to check that the mass-assembly his- 
tory (for fixed halo mass) does not shift the energy-infall relation 
much, we suspect that this will not be the major driver for inflating 
the width of the infall-time PDFs. This is because the Milky Way is 
likely to have had a relatively quiescent merger history for the ma- 
jority of its life. Simulations of disk galaxies have indicated that the 
thinness of the Milky Way stellar disk cannot be reproduced if the 
Milky Way had a cosmologically- typical 10:1 merger since 2 = 1 
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dPurcell. Kazantzidis & Bullockl [20091) . As we showed in Sec. [3] 
the gravitational potential in the outskirts of halos does not change 
much with time if halos are in the "slow growth" phase charac- 
teristic of quiescent halo evolution. Since the halo potential does 
depend on the overall mass distribution of the dark-matter halo, we 
expect the halo mass rather than the mass-assembly history to drive 
the normalization and the slope of the energy-infall relation. 

The main reason to simulate many Milky Way-type galaxies, 
even if there were no uncertainty in the Milky Way gravitational 
potential and its evolution, is to get better statistics on the subhalos. 
Since VL2 has a relatively small number of halos, we had to choose 
relatively large windows in r and V r to select subhalo samples for 
each Milky Way dwarf, and we only had samples of 20 subhalos for 
each dwarf (Sec. |4j. If we had a much larger subhalo sample from 
which to work, we could select only subhalos with r and V r within 
the observational uncertainties and have a large enough sample to 
create a smooth infall-time PDF. 

On the observational side, the main way to improve the infall- 
time PDFs is to include proper-motion data. As we have seen in 
Sec. [4] proper-motion measurements sharpen the PDFs, but only 
if the uncertainty in the proper motion translates to ~ 50 km s -1 . 
However, in the case of Fornax, there is still a double-peaked PDF 
in spite of an excellent proper motion measurement. It is possible 
that better subhalo statistics will sharpen the PDF since the win- 
dow in r and V r from which we select a representative simulated 
subhalo sample for each dwarf will narrow. 

The best thing to do, in terms of producing robust infall-time 
PDFs, is to simulate a number of Milky Way-type systems with 
different halo masses and with baryons. This will give us both good 
subhalo statistics and an idea of how much the dark and baryonic 
masses affect the energy-infall relation since there are significant 
uncertainties in both for the Milky Way. It would also be useful to 
get or improve proper-motion measurements of the observed dwarf 
galaxies. 

5.3 Other consequences: stellar streams and dark matter 
flows 

We have focused on the energy-infall relation for surviving dwarf 
galaxies, but we may also apply this relation to their tidal debris and 
to those dwarfs that did not survive. In particular, we expect there 
to be an energy-infall relation for tidal streams. Unless dynamical 
friction or satellite-satellite encounters are important for a satellite, 
its energy should not change much after infall onto the Milky Way 
(as we showed in Sec.fJJ. This means that tidally stripped material 
should "remember" the satellite's infall energy. 

Since the shape of the tidal stream approximately marks out 
the orbit that the surviving part of the satellite should have had, 
one can estimate the energy of the progenitor satellite based on 
the shape of the tidal stream. The next step would be to find infall- 
time PDFs for the progenitor satellite and compare that to the abun- 
dance pattern and inferred star-formation history of the stars in the 
streams. Though more work needs to be done in order to check this 
possibility, it could prove to by particularly important for the ef- 
forts of near-field cosmology. Specifically, the progenitor galaxies 
that formed the earliest, at times approaching the epoch of reoiniza- 
tion, are more likely to be accr eted early and be shredded by now 
(e.g. iBullock & Johnstonll2005l) . If one of the goals of Milky Way 
studies is to use local observations to constrain high-redshift phe- 
nomena, it would be useful to be able to associate stellar streams 
with accretion (and ultimately, formation) times in as many ways 
as possible. 



The energy-infall relation also explains t he typical energy 
of the dwarf galaxy debris field discussed by llisanti & Spergell 
J201 1). Those authors are concerned with characterizing the dark- 
matter particle population at the Sun's position in the Milky Way 
from the tidal debris of the subhalos surviving today that had identi- 
fiable progenitors in VL2 at z — 9, which they define as the epoch 
of reionization. They find that the this debris has a characteristic 
velocity at the Sun's position. We would explain this characteristic 
velocity feature as arising from the energy-infall relation. Since the 
Sun is deep within the Milky Way's potential well, we find from 
Fig. [3] that any surviving subhalos (and their associated tidal de- 
bris) should have been accreted at early times. In other words, the 
infall-time PDF of surviving subhalos and their streams that might 
pass through the Sun's orbit is sharply peaked at high tinfall- This 
means that the typical orbital energy today £ , and hence the typical 
dark-matter particle speed, is also peaked. 



6 SUMMARY AND CONCLUSIONS 

We used a single simulatio n of a Milky- Way-like ha l o, the 
Via Lactea II (VL2) halo foiemand, Kuhlen & Madaul 120071 : 
iDiemand et alJl2008l ; lKuhienll2010l) . to investigate the relationship 
between the kinematic properties of subhalos and their infall times. 
We found that there exists a tight correlation between the binding 
energy of a subhalo to the host at z = and the time at which the 
subhalo last passed inward through the host virial radius (see Fig. 
1). We found that the origin of this "energy-infall relation" was the 
on-average near conservation of the energy of the subhalo from the 
time of infall. This means that the orbits of the subhalos are deter- 
mined by the mass-assembly history of the host, since the energy 
of the infalling subhalo must be similar to that of a particle on a 
circular orbit at the virial radius at the time of infall. 

Assuming that the Milky Way has a similar mass-assembly 
history and gravitational potential as the VL2 halo, we were able to 
assign infall-time probability distribution functions for the known 
Milky Way satellite galaxies based on their present-day kinematics 
and the energy-infall relation. We found reasonably peaked infall- 
time PDFs even in galaxies where proper motions were not avail- 
able (and hence only had an upper limit on the binding energy from 
their radii and radial velocities), but the addition of proper motion 
constraints sharpen the PDFs considerably. For example, the or- 
bital energies of Carina, Ursa Minor, and Sculptor are all strongly 
indicative of early accretion, more that 8 billion years ago. Con- 
versely, Leo T, Leo I, and the LMC were all recently accreted, 
within the last few billion years. Fornax, Leo II, and the ultrafaint 
dwarf Canes Venatici I are all examples of intermediate lookback 
time accretions. Accretion time estimates for each dwarf are pro- 
vided in Table 1 . 

When comparing the infall-time PDFs for individual Milky 
Way satellite galaxies with their inferred star-formation histories, 
we found a number of interesting trends. For the classical dwarf 
galaxies, we found that the infall time occurred at a similar time as 
the end of star formation in about half the sample, but that star for- 
mation continued for a short while before ending in the other half. 
Several of the dwarfs had at least one burst of star formation imme- 
diately after infall. For the ultrafaint dwarf galaxies, we found that 
about half the galaxies fell in quite early and could not determine 
whether star-formation quenching occurred before or after infall. 
Several of the dwarfs showed some evidence for a small amount 
of star formation after infall (Leo IV, Canis Venatici L and possi- 
bly Ursa Major II), but Hercules clearly stopped forming stars long 
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before it fell into the Milky Way. In those cases, the bulk of star 
formation was quenched long before the dwarfs became satellites 
of the Milky Way. The diversity of the offsets, both in magnitude 
and sign, in the infall time versus the end of star formation sug- 
gests that quenching is not a single and uniform process for these 
smallest of galaxies. It means that any star-formation-density or 
morphology-density relation for the dwarf galaxies does not share 
the same origin as those relations for more massive galaxies. 

The energy-infall relation that we have presented here will 
hopefully inspire explorations that benefit from a larger number of 
simulations, including eventually baryonic physics. If the correla- 
tion between energy and infall time proves to be universal, even in 
some renormalized way, it will provide a new avenue for testing for- 
mation scenarios for satellite galaxies. In this work we have taken 
the first tentative steps in that direction, part of the larger goal to 
put these interesting objects to use as the Rosetta Stones of galaxy 
formation that they have longed promised to be. 
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